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Abstract. The goal of the DIRAC experiment at CERN is to measure with high 
precision the lifetime of the 7r + 7r~ atom (A 27r ), which is of order 3x 10~ 15 s, and thus to 
determine the s-wave 7T7r-scattering lengths difference |oo — 0,2]- A-i-n atoms are detected 
through the characteristic features of 7r + 7r~ pairs from the atom break-up (ionization) 
in the target. We report on a first high statistics atomic data sample obtained from 
p Ni interactions at 24 GeV/c proton momentum and present the methods to separate 
the signal from the background. 
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Dedicated to the memory of Lucien Montanet 

1. Introduction 

The measurement of the lifetime of the atom A 2n [1], which is essentially 

determined by the 7t + tt~ — > 7r°7r° reaction, enables the determination of the combination 
\a — a 2 \ of the s-wave 7T7r-scattering lengths for isopins 1 = and 2 in a model 
independent way 013111130111111 according to [7] 
1 2 

— = - a 3 p \a - a 2 \ 2 (1 + 5), (1) 

Tn y 

with Tis the lifetime of the atomic ground state, a the electromagnetic coupling constant, 
p the 7T° momentum in the atomic c.m.s., and 5 = (5.8 ± 1.2) x 10~ 2 a correction due 
to QED and QCD 0. 

The 7T7T scattering lengths Oq and a 2 have been calculated within the framework of 
Chiral Perturbation Theory |Hj by means of an effective Lagrangian with a precision 
at the percent level [TU]. The lifetime of A 2n in the ground state is predicted to 
be TU = (2.9 ±0.1) x 1CT 15 s. These results are based on the assumption that 
the spontaneous chiral symmetry breaking is due to a strong quark condensate as 
recently suggested [TTjlT2"] . An alternative scenario with an arbitrary value of the quark 
condensate ^3] allows for larger do, a 2 compared with those of the standard scheme |l()j . 
A measurement of the scattering lengths will thus contribute crucially to the current 
understanding of chiral symmetry breaking in QCD and constrain the magnitude of the 
quark condensate. 

The differential production cross section of A 2n atoms can be obtained from the 
double differential two-pion production cross section (cjo) without a Coulomb final state 
interaction 



llfT ' } (2vr) 3 ^ |* n (0)| 2 da ° 



(2) 

Pl=P2 



dp a Ma dpidp 

where the differential production cross section <r„ for atoms with principal quantum 
number n and zero angular momentum depends on mass (Ma), momentum (pa) and 
energy (Ea) of the atom in the lab frame, and on the square of the Coulomb atomic 
wave function at zero distance |\l/ ra (0)| 2 . The laboratory momenta of the 7r + and ir~ are 
denoted by pi, p 2 , respectively. On the basis of Eq.EJand using the Fritiof 6 generator, 
yields for A 2n in proton nucleus interactions have been calculated as a function of their 
energy and angle in the proton energy range from 24 GeV to 1000 GeV [THll5| ITo]. For 
a Ni target and a 24GeV/c proton beam, ~ 7 x 10~ 7 atoms are produced per proton 
interaction, of which, however, only ~ 1 x 10~ 9 are detectable in the experiment, due 
to momentum and angular acceptance of the DIRAC apparatus. 

The method for observing A 2n and measuring its lifetime has been proposed in |14j . 
Pairs of 7r + 7r~ are produced mainly as unbound ("free") pairs, but sometimes also as 
A 2n . The latter may either decay, get excited or break up into 7r + 7r~ pairs (atomic pairs) 
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after interacting with target atoms. Due to their specific kinematical features, these 
atomic pairs are experimentally observable. For thin targets (10 _3 X ) the observable 
relative momentum Q in the atomic pair system is Q < 3MeV/c. Their yield is 
~ 10% 4- 20% of the number of free pairs in the same Q interval. In Fig.^the relative 
momentum distributions of atomic pairs are shown for Q and Ql (relative momentum 
component along the flight direction of the pair) at the moment of break-up and at 
the exit of the target, as obtained with an event generator based on J7j. The number 
of atomic pairs is a function of the atom momentum and depends on the dynamics of 
the A 2n interaction with the target atoms and on the A 27T lifetime ^H]- The theory 
of the A 2n interaction with ordinary atoms allows the calculate of the relevant cross 
sections |T7l HH HH1 1201 |2H 1221 12H I2SI I2ZI I2HJ - For a given target thickness, the 
theoretical breakup probability for A 2n is precise at the 1% level and uniquely linked to 
the lifetime of the atom JX^EHJ - In Fig. E] the break-up probability as a function of the 
lifetime is displayed for a 94 /im thick Ni target and for atomic pairs accepted in the 
DIRAC spectrometer. 
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Figure 1. Relative momentum 
distributions {Q, Ql) for atomic 
7r + 7r~ pairs at break-up and at the 
exit of the target (event genera- 
tor, no spectrometer simulation). 
Note that Ql is almost not af- 
fected by multiple scattering in 
the target. 



Figure 2. Break-up probability 
Pbr as a function of the ground 
state lifetime t\ s of the Am 
atom for a Ni-target and DIRAC 
conditions. 



The first observation of A 2n [HO] has been achieved in the interaction of 70 GeV/c 
protons with Tantalum at the Serpukhov U-70 synchrotron. In that experiment, the 
atoms were produced in an 8 /im thick Ta target, inserted into the internal proton beam. 
With only 270 ± 50 observed atomic pairs, it was already possible to set a lower limit 
on the A 2n lifetime jSBE2]: r > 1.8 x 10~ 15 s (90% CL). 



Detection of 7v + tt atoms at DIRAC 



5 



In this paper, we present the first high statistics experimental data on A^^ 
production on a Ni target at an external proton beam of the CERN PS and demonstrate 
the feasibility of the lifetime measurement. We are not attempting to deduce a lifetime 
here, as this requires a highly involved analysis of the normalization when evaluating 
the break-up probability. 

2. The DIRAC experimental setup 

The DIRAC setup is designed to detect oppositely charged pion pairs of low relative 
cm. momenta with high resolution using a magnetic double arm spectrometer at a 
24 GeV/c extracted proton beam of the CERN PS and an especially low material budget 
in the secondary particle path. The spectrometer setup is shown in Fig.|31 A detailed 
description may be found in [33] . 




Mil 



Figure 3. Schematic top view of the DIRAC spectrometer. Upstream of the magnet: 
target, microstrip gas chambers (MSGC), scintillating fiber detectors (SFD), ionization 
hodoscopes (IH) and iron shielding. Downstream of the magnet: drift chambers (DC), 
vertical and horizontal scintillation hodoscopes (VH, HH), gas Cherenkov counter (Ch), 
pre-shower detectors (PSh) and, behind the iron absorber, muon detectors (Mu). 

The proton beam intensity during data taking was 0.9 x 10 11 per spill with a 
spill duration of 400 -r- 450 ms. The beam line was designed such as to keep the 
beam halo negligible [33]. The horizontal and vertical widths of the beam spot were 
o"hor = (0.80 ± 0.08) mm and <7 vert = (1.60 ± 0.07) mm, respectively [SSI- The targets 
we report here were 94 and 98 /im thick Ni foils, corresponding to ~ 6.4 x 10~ 4 nuclear 
interaction probability [36J or 6.7 x 10~ 3 radiation length. The transverse dimensions of 
the circular targets (4.4 cm diameter) were sufficient to contain the proton beam fully 
and to exclude possible interactions of beam halo with the target holder. 

The spectrometer axis is inclined upward by 5.7° with respect to the proton beam. 
Particles produced in the target propagate in vacuum up to the upstream (with respect 
to the magnet) detectors and then enter into a vacuum chamber which ends at the 
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exit of the magnet. The exit window is made of a 0.68 mm thick Al foil. The angular 
aperture is defined by a square collimator and is 1.2 msr (±1° in both directions). The 
dipole magnet properties are B = 1.65 T and BL = 2.2 Tm. 

The upstream detectors are 2.5 to 3 m away from the target and cover an area of 
roughly 10 x 10 cm 2 . The microstrip gas chambers (MSGC) consist of four planes: X, Y, 
U and V, with rotation angles of 0°, 90°, 5°, 85° with respect to the X-plane. Each plane 
has 512 anode strips with a pitch of 0.2 mm. Clustering results in a spatial resolution for 
single tracks of 54 /im. The two planes (X, Y) of the scintillating fiber detector (SFD) 
provide both coordinate and timing information. Each plane contains 240 fiber columns 
(0.44mm pitch), each column consisting of 5 fibers of 0.5 mm diameter. They are 
read out through multichannel position-sensitive photomultipliers and an analog signal 
processor that produces information for the appropriate TDC channel. The space and 
time resolutions are 130 /im (rms) and 0.8 ns, respectively. The analog signal processor 
merges two adjacent hits into one, depending on the relative pulse hight of the two signals 
and on their time difference. No merging takes place if the time difference is larger 
than 5 ns. The ionization hodoscopes (IH) serve for fast triggering and identification 
of unresolved double tracks through an energy loss measurement. The IH detector is 
composed of two vertical (X) and two horizontal (Y) layers, each with 16 slabs of plastic 
scintillator (1 x 7 x 110mm 3 ). The read-out provides logic (TDCs, trigger processors) 
and analog (ADCs) information. The total thickness of all upstream detectors, including 
the vacuum channel windows, is 6.0 x 10~ 2 X . 

The two arms of the spectrometer are identically equipped. Four sets of drift 
chambers are used (DC1 to DC4). DC1 and DC4 have two X and two Y planes each, 
DC1 has in addition two inclined W planes (11.3° with respect to the X-wires). DC2 
and DC3 have one X and one Y plane each. The sensitive areas range from 0.8 x 0.4 m 2 
to 1.28 x 0.4 m 2 , the signal wire pitch is 10 mm. The space resolution is better than 
90 /im. The vertical hodoscopes (VH) supply time-of-flight information and serve trigger 
purposes. The hodoscope is made of 18 plastic scintillation counters (2.2 x 7 x 40cm 3 ). 
The time resolution is 127 ps. The horizontal hodoscope (HH) consists of 16 plastic 
scintillators (2.5 x 2.5 x 130 cm 3 ). It serves essentially trigger purposes. The threshold 
Cherenkov detectors (Ch), are used to identify electrons (positrons) and to reject pairs 
containing an electron and/or a positron. The radiator is Nitrogen gas at normal 
pressure and ambient temperature. The average number of photoelectrons for particles 
with (3 ~1 is larger than 16 and the efficiency more than 99.8%. Pion contamination 
above the detection threshold is estimated to be less than 1.5%. The pre-shower detector 
(Psh) consists of 8 elements, each comprising a Pb converter and a scintillator. Off-line 
analysis of the amplitudes from Psh provides additional e/V separation. Each muon 
detector is made of a thick iron absorber followed by two planes of plastic-scintillator 
counters with 28 counters per plane. 

The momentum range covered by the spectrometer is 1.2 4- 8 GeV/c. The relative 
resolution on the lab-momentum is dominated by the multiple scattering downstream 
of the vertical SFD detector and the spatial resolution of the drift chambers. While 
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the former leads to a momentum independent resolution, the second leads to a slight 
increase with momentum. Studies of the A signal in 7r~p track pairs [37j result in 
o~ p /p = 2.8 x 10~ 3 and increasing to 3.3 x 10~ 3 at 8 GeV/ c. The resolutions on the relative 
c.m.-momentum of 7r + 7r~ pairs from atomic break-up, Q = (Q x , Q y , Ql), are in the plane 
transverse to the total momentum pJ+%-, ctq x « oq v < 0.49MeV/c and longitudinally 
°Ql = 0-50 MeV/c fE] t- The setup allows to identify electrons (positrons), protons with 
p < 4 GeV/c (by time-of-flight) and muons (cf. Fig. 23 in [33 ). It cannot distinguish 
7r- from if-mesons, but kaons constitute a negligible contamination [HUES!- 

The trigger system [40J reduces the event rate down to a level acceptable for 
the data acquisition system. The on-line event selection keeps almost all events with 
Ql < 22MeV/c, Q T = Jof+Q^ < lOMeV/c and rejects events with Q > 15MeV/c 
progressively |4*T] . In the first level trigger [12] a coincidence of VH, HH and PSh signals 
and anticoincidence with Ch signals in both arms is treated as a pion-pair event. A 
condition on the HH of the positive and negative arms, |HH + — HH_| < 2 slabs, rejects 
events with Q y > 10 MeV/c. Fast hardware processors [4THl43p44jl43] are used to decrease 
the first level trigger rate by a factor 5.5. The trigger rate in standard conditions was 
around 700 per spill. The trigger accepts events in a time window ±20 ns with respect 
to the positive arm and thus allows for collecting also accidental events. The trigger 
system provides parallel accumulation of events from several other processes needed for 
calibration such as e + e~ pairs, A — > pn^ or 7r ± 7r + 7r _ final states. 

The data acquisition jlHlllZI accepted up to 2000 events per spill, at spill intervals 
as short as one second. 

Dead times were studied in a run with 30% higher intensity than normal and depend 
on average beam intensity as well as on micro duty cycle of the beam. An overall dead 
time of 33% was found for triggers and data acquisition, with an additional 10% due 
to front-end electronics [IB]- At normal intensities, dead times are lower. Biases due 
to dead times could not be found. When selecting data offline for further analysis, run 
periods, have been eliminated where problems with detectors or spill structure or micro 
duty cycle were found. 

The full setup, including detectors, detector responses, read-outs, triggers and the 
magnet has been implemented into the detector simulation DIRAC-GEANT [49J and 
into the DIRAC analysis package ARIANE (HH] such that simulated data can be treated 
in the same way as real measurements. 

3. Track reconstruction 

Events of interest consist of two particles of opposite charge with very low Q, resulting 
in two close lying tracks upstream that are separated downstream into the two arms of 
the spectrometer. 

The magnetic field was carefully mapped and fine-tuned to a relative precision 
of 10~ 4 . The final field map is used in the Monte Carlo GEANT simulation of the 

| A momentum dependence of ~ 0.016 [MeV/c] x p + _ (p in GeV/c) is included in these numbers. 
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experiment |1H]. For tracking, the field is summarized by a transfer function that uses 
as input the position, direction and momentum of a track in a reference plane at the exit 
of the magnet and calculates the position and direction of that track in a reference plane 
five meters upstream of the magnet. The transfer function consists of 4 polynomials with 
5 parameters each. 

Track reconstruction starts from the downstream part. A track candidate is 
searched for in the horizontal and vertical planes of the drift chambers separately. A 
"horizontal" candidate must have a hit wire in one of the horizontal planes of the first 
(DC1) and last (DC4) chamber, as well as corresponding hits in the HH-detectors. They 
define a straight line. Hits close to that line are looked for in the other four horizontal 
planes of DC1 to DC4. Based on an acceptance window, a candidate is accepted if at 
least two more hits are found. The "vertical" candidates are found analogously. Finally, 
candidates from horizontal and vertical wires are matched using the inclined planes. 
Projecting the spatial track parameters by means of the magnetic field polynomials 
to the center of the beam spot at the target provides a rough momentum estimate 
for each track candidate. Using this estimate and the precise position measurements 
deduced from drift-times, an overall track-fit for each candidate is made using a standard 
least-squares method, where the full error matrix is used, including the correlations (off- 
diagonal terms) induced by multiple scattering. Candidates are retained on the basis 
of a confidence level cut. The downstream track candidates are thus defined spatially 
with high precision. This is demonstrated in Fig.HJ where the difference between the 
measured hit coordinate in the first DC1 X-plane and the coordinate obtained with the 
track parameters at the same DC1 plane is shown. The spatial and angular resolutions 
depend on the particle momentum. They show a 1/p dependence and level off at 
6GeV/c. As an example the angular resolution in X direction is shown in Fig. 03 The 
one in the Y direction is about 4% better. 

The multiple scattering in the upstream detectors is such that in practice they have 
to be considered as the effective source of each track, and not the target as supposed 
above. Therefore, each of the above tracks, projected through the magnet onto the 
beam spot defines a reference line. Along this line, hit candidates are searched for in 
the upstream detectors within spatial windows defined mainly by multiple scattering 
and within time windows defined by the time-of-flight from the upstream detectors to 
the vertical hodoscopes downstream and their time resolutions. 

Using all identified upstream information, a track-fit is made for each track- 
candidate by means of the Kalman-filter method [51] . starting from the first downstream 
hit and ending at the exit window of the vacuum tube upstream of the MSGC detector. 

With the hypothesis that the track originates from the target, the intersection of 
the proton beam with the target provides another measurement point for the Kalman 
filter, whose uncertainty is given by the measured intensity distribution of the proton 
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Figure 4. Difference between 
a hit coordinate measured by the 
first X-plane of DC in the right 
arm and a coordinate calculated 
with track parameters. 



Figure 5. Resolution for angles 
of DC-tracks in X projections 



beam across the target (cf. section EJ) §. A cut on the distance (15 mm) of the track 
from the beam spot leads to the final track selection 

Alternatively, assuming that a track pair originates from the same interaction, each 
pair candidate is fitted with the constraint that both tracks intersect in the central plane 
of the target ("vertex fit"). The complete 5x5 uncertainty matrices of the two tracks 
are used. Tracks originating far from each other are rejected by a threshold confidence 
level. This procedure provides track parameters which are independent of the precise 
knowledge of the beam position and beam width. For the subsequent analysis these two 
procedures do not use the MSGCs, for reasons of optimum efficiency 

Full tracking uses MSGC and SFD to connect all upstream hits by straight lines. 
These track candidates are then matched with the downstream candidates. A vertex 
confidence level is calculated for each track pair, corresponding to the hypothesis that 
both tracks intersect at a common point lying on the target plane. Correlations from 
the estimated errors from multiple scattering are fully taken into account [53]. Both 
tracks are refitted at the end, under the constraint of having a common vertex [S^j . 

To illustrate that the reconstructed events do in fact originate from the target, the 
distances between two tracks at the target are shown in Fig-EJ No vertex fit was done 
for these data. 

The correctness of the alignment and magnetic field description were verified 
by studying the decay of A particles jHH]- Events with one downstream track per 
spectrometer arm and a time-of-flight difference between the positive particle and the 
negative one between and 1.3 ns were selected. Fig.[7| shows, for a typical control 

§ The beam position is continuously monitored and calibrated during data taking. 



Detection of tt + tt atoms at DIRAC 
12000 H 



10 



10000 



8000 



6000- 



4000- 



2000 



RMS_x= 3.24 mm 
RMS_y= 3.64 mm 

solid line: x 
dashed line: y 



-2-1012 

Ax,y [cm] 

Figure 6. Distances Ax and Ay 
between two reconstructed tracks 
at the target in the transverse 
plane. 



M A = 1115.790(31) MeV/ c 2 
o M = 0.395(07) MeV/c2 




1112 1116 



1120 1124 1128 

M pK - [MeV/ C 2] 



Figure 7. Invariant proton-pion 
mass (M p7T -) distribution. The 
dashed line represents a fit of a 
Gaussian plus a straight line. 



data sample, the invariant mass distribution from the two tracks under the hypothesis 
proton-pion. A fit of a Gaussian on a linear background yields Ma = (1115.790 ± 
0.031 siai ) MeV/c 2 with a MA = (0.395±0.007 siai ) MeV/c 2 . The measured width is entirely 
due to track reconstruction. The difference (M p7r - — M^PDG jSHl)) in units of the 
reconstruction error provides a measure of the correctness of the error estimation. The 
resulting distribution is fitted by a Gaussian with a = 1.028 ± 0.006, showing that the 
reconstruction errors are slightly underestimated. The same distribution was studied as 
a function of the pion momentum and found to be independent of it. The long term 
stability of the apparatus has been controlled using the A mass and the corresponding 
widths. 

4. Selection criteria 

An event is rejected if more than two downstream tracks in either of the two arms are 
reconstructed. For each track the associated VH and HH hodoscope slabs are required. 
In the case of two tracks in one arm, the earlier in time is taken for further analysis. 
Events with more than one track per arm constitute less than 4% of the event sample. 

Horizontal and vertical SFD hit candidates to be associated with a track (see 
section El for first momentum estimation) must have times within a window of ±4 ns 
(corresponding to ±5<t) with respect to the associated VH slab. Moreover, the 
hit candidates must be found in a spatial window of ±(0.2 + 4.8/p [GeV/c]) cm 
(corresponding to ±6<r) with respect to the point of intersection of the track candidate 
with the SFD plane. This window is defined by multiple scattering in the downstream 
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material when backward-extrapolating the track. At least one hit candidate is required. 
If there are more than four, the four closest to the window center are kept. 

After the first stage of Kalman filtering, only events are kept with at least one 
track candidate per arm with a confidence level better than 1% and a distance to the 
beam spot in the target smaller than 1.5 cm in x and y. For such tracks, the SFD hits 
are redetermined in a narrower (±lcm) window around the new track intersect. Events 
with less than four hits in the window per SFD plane and not more than six in both 
planes are retained. Finally, all track pairs with \Q X \ < 6MeV/c, \Q y \ < 6MeV/c and 
\Ql\ < 45MeV/c are selected for further analysis. 

This preselection procedure retains 6.2% of the titt data obtained with the Ni target 
in 2001. The time difference between the positive and the negative arm, measured by 
the VHs, for events that passed the preselection criteria is shown in Fig.lHl 

For the final analysis further cuts and conditions are applied: 

• "prompt" events are defined by a time difference (corrected for the flight path 
assuming pions) measured by the VHs between the positive and the negative arm, 
|At| < 0.5 ns, corresponding to a ±2.7<r cut. 

• "accidental" events are defined by time intervals —15 ns< At < —5 ns and 
7ns< At < 17 ns, conditioned by the read-out features of the SFD detector (cf. 
section^ and suppressing protons (cf. Fig.[HJ). 

• protons in prompt events are rejected by requiring momenta of the positive particle 
to be P+ < 4 GeV/c. 

• e ± and are rejected through appropriate cuts on the Cherenkovs, the Pre-shower 
and the Muon counters [5T>] . 

• Q T < 4MeV/c and \Q L \ < 22MeV/c. The Q T cut preserves 98% of the atomic 
signal, the Ql cut preserves background outside the signal region for defining the 
background below the signal. 

• the vertex-fit with highest confidence level from track and vertex fits is retained. 

• only events with at most two preselected hits per SFD plane are accepted. This 
provides the cleanest possible event pattern. This criterion is not used with the full 
tracking (cf. section HH). 

5. Signal and background in 7r + 7r atom detection 

Once produced in proton-nucleus interactions, all atoms will annihilate if they 
propagate in vacuum. In a finite target, however, they interact electromagnetically with 
the target atoms, and some of them break up. The event generator for 7r + 7r~ pairs 
from A 27T break-up yielded Q and Ql distributions shown in Fig.[T] (cf. section [IJ. At 
the exit of the target, multiple scattering has caused a considerable broadening of the 
Q-distribution while the Ql distribution has remained unchanged. The Ql distribution 
at break-up is already much narrower than the Q-distribution because the break-up 
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Figure 8. Time difference be- 
tween positive and negative ver- 
tical hodoscope slabs associated 
with the selected tracks. The 
asymmetry for positive differences 
is due to time correlated protons. 



Figure 9. Difference in hit 
number for selected SFD hits in 
the X-plane, for accidental and 
prompt events. Note the effect 
of merging of adjacent hits for 
prompt events. Arbitrary scale. 



mechanism itself affects mostly Qt (analogous to 5 electron production). These features 
led us to consider in the analysis not only Q but also Ql distributions . 

The atomic pairs are accompanied by a large background. The Q distribution of 
all 7r + 7r~ pairs produced in single proton-nucleus interactions in the target is described 
by: 

dN/dQ = dN c /dQ + dN nC /dQ + dn A /dQ. (3) 

Nq is the number of 7r + 7r _ pairs originating from short-lived sources (fragmentation, 
rescattering, mesons and excited baryons that decay strongly). These pairs undergo 
Coulomb interaction in the final state (Coulomb pairs). N n c is the number of 7r + 7r~ 
pairs with at least one particle originating from long-lived sources (mesons and baryons 
that decay electromagnetically or weakly). They do not exhibit Coulomb final state 
interaction (non-Coulomb pairs). Finally, ua is the number of 7r + 7r~ pairs from Ai-k 
break-up. 

Accidental 7r + 7r~ pairs (dN acc /dQ) originate from different proton-nucleus 
interactions and are uncorrelated in time, i.e. they are neither affected by Coulomb 
nor by strong interaction in the final state. Such events may also belong to the time 
window that defines prompt events (see section HJ). 

These backgrounds, as they are measured and reconstructed, are needed for 
subtraction from the measured data in order to obtain the excess produced by the 
atomic signal. The backgrounds may be obtained in different ways. 

One method is based on a Monte Carlo modelling of the background shapes 
using special generators for the non-Coulomb (nC) and Coulomb (C) backgrounds 
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and accidental pairs Uniformity in phase space is assumed for the backgrounds 

(dN/dQ ~ Q 2 )-, modified by the theoretical Coulomb correlation function A i Q eo {Q) [SH] 
for C-background. Effects of the strong final state interaction and the finite size of the 
pion production region [59 on the Coulomb correlation function are small and neglected 
here. The differences of pion momentum distributions for the different origins of 7r + 7r~ 
pairs are implemented (cf. (HOIEZI ) - The generated events are propagated through the 
detector (GEANT [H5J). Simulating the detectors, detector read-outs and triggers and 
analyzing the events using the ARIANE analysis package of DIRAC [SD] result in the 
high statistics distributions dN c /dQ, dN nC /dQ (cf. Eq.^J) , and dN &cc /dQ. These 
distributions are used later to analyze the measured distributions. Though not necessary 
for the signal extraction, the A 27T signal was also simulated yielding dnj^jdQ. Monte 
Carlo simulated distributions were generated with about 10 times the statistics of the 
measured data. The simulated backgrounds were obtained with dedicated generators, 
without any additional tracks from the proton nucleus interaction. A special Monte 
Carlo simulation with additional background tracks leads essentially to a reduction in 
reconstruction efficiency. 

The second method circumvents the Monte Carlo simulations of background, 
detectors and triggers by relating the backgrounds to the measured accidental 
background. It was successfully applied in (30] . The assumption is that accidentals can 
be used to describe the distribution of free 7r + 7r~ pairs, which then must be corrected for 
final state interactions. By denoting $(Q) = dN^ as /dQ, the experimental correlation 
function R{Q) is given by: 



Here Ac(Q) is the Coulomb enhancement function j?8- smeared by multiple 
scattering in the target and the finite setup resolution. In practice it can not be 
obtained by folding procedures but must be obtained from Monte Carlo simulations 
described above. N and / are free parameters, and $(Q) is the measured spectrum 
of accidentals corrected for the differences between prompt and accidental momentum 
distributions 60 . 

The measured accidental distributions had to be corrected for the different recording 
conditions as compared to prompt events, such as the SFD merging of adjacent hits, or 
track identification by time, which is impossible for prompt events in case of ambiguities. 
In Fig.lHlthe difference in selected SFD hits is shown for accidental and prompt events for 
the X-plane (the Y-plane is similar). The merging feature of the SFD read-out together 
with a ~ 5% single track inefficiency of the SFD lead to the enhancement for A-SFD=0, 
and to the dips left and right of it. The correction of the measured accidental spectra for 
prompt conditions was done in two ways. On an event by event basis adjacent hits of the 
measured accidentals were merged into one hit according to the measured probability, 
and that hit was given the times of the two original hits. Then tracking was started, 
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and the timing conditions could be applied as described above (cf. section HJ. An other 
way of correcting was to construct accidental distributions from uncorrected events, but 
giving up the timing conditions. The resulting distributions were then given weights 
according to the measured merging probabilities. 



6. Experimental data and atomic signal extraction 

DIRAC began data taking in autumn 1999. Here we present data taken in 2001 using 
94 /xm and 98 /im thick Nickel targets. The integrated proton flux through the target 
for these data is 8.6 x 10 16 , corresponding to 5.5 x 10 13 pNi interactions and to 6.4 x 10 8 
recorded irn triggers. The reconstructed accidental and prompt tttt pairs (see Section 0] 
for the cuts) are shown in Fig.fTUl 
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Figure 10. Experimental Q and Ql distributions for prompt and accidental rnr events. 
Continuous lines serve as guides for the eye. The sharp peak at Q = 4MeV/c is due 
to the cut Qt < 4MeV/c. The number of events is given for the displayed momentum 
window, the bin- width is 0.5MeV/c. Note the strong Coulomb enhancement in the 
Ql distribution for prompt events. 



The Monte Carlo backgrounds or the correlation function R(Q) (cf. Section EJ) are 
fitted to the measured prompt spectra in Q and Ql intervals which exclude the atomic 
signal (typically Q > 4MeV/c, Ql > 2MeV/c). The amount of time-uncorrelated 
events in the prompt region was determined to be 6.5% of all prompt events by 
extrapolating the accidental At distribution to zero and was subtracted from the prompt 
distributions (cf. Fig.IHJ) • The fits provide the relative amounts of non-Coulomb and 
Coulomb backgrounds and the parameters N and / of Eq.0] As a constraint, these 
background components must be the same for Q and Ql- 

The analysis using Monte Carlo backgrounds and no vertex fit is summarized in 
Fig.fTIl The background composition is indicated, and the excess at low Q and Ql is 
clearly seen. Subtraction provides the residuals, also shown in Fig. El which represent 
the atomic signals. As expected from Fig.^ the signal is narrower in Ql than in Q. 



Detection of tt + tt atoms at DIRAC 



15 



The simulated shape is in agreement with the data. The signal strength has to be the 
same in Q and Ql if the background is properly reconstructed. The observed difference 
demonstrates that the backgrounds are consistent at the per mille level. The fact that 
outside of the signal region the residuals are perfectly zero demonstrates the correctness 
of the Monte Carlo simulation. 
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Figure 11. Top: Experimental Q and Ql distributions after subtraction of the 
time uncorrelated background, approximated with Monte Carlo backgrounds (dashed 
lines). Bottom: Residuals after background subtraction. The dashed lines represent 
the expected atomic signal shape. Notice that the signal strengths in Q and Ql are 
about the same. The bin-width is 0.25MeV/c. 



The results of the analysis using the accidentals as a basis for background modelling 
are shown in Fig.^J The clear deviation of R(Q) at low Q from unity is due to the 
attractive Coulomb interaction in the final state. There is good agreement between 
the experimental and fitted "correlation" function R for Q > 2.5MeV/c, whereas for 
Q < 2.5MeV/c the experimental correlation function (left hand side of Eq.0J shows 
significantly higher values, due to the presence of atomic pairs. In order to extract the 
number of atomic pairs the fit function (right hand side of Eq.0J) has been multiplied by 
the accidental distribution $(C}). The subtraction of the latter from the full measured 
spectrum leads to the atomic pair signal in Fig.lT2l 

Based on the same data sample of preselected events, the analysis using Monte 
Carlo backgrounds and no vertex fit was repeated allowing for three hits instead of only 
two in the search window (cf. section 0} and selecting the hit closest to the search 
reference trajectory. This yields about 25% more events in the signal, but 43% more in 
total (see Tabled)), thus the signal/background quality diminishes. 

An analysis fully independent on the upstream detectors was done using only 
information from the downstream detectors. Events were selected from raw data along 
with a looser cut in Qt- The particles were assumed to be emitted from the center of 
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Figure 12. Top: Experimental Q and Ql correlation functions (experimental points). 
The background correlation function is shown as a smooth curve. Bottom: Residuals 
after background subtraction. The dashed lines represent the expected atomic signal 
shape. The bin- width is 0.5MeV/c. 



the beam spot at the target. For this analysis only the Ql distribution shows a signal, 
because multiple scattering in the upstream detectors does not distort Ql so much but 
destroys Qt and hence Q. With respect to the standard selection the signal increases 
by about 40-50%. The total number of events is, however, almost 3.5 times larger (see 
Table^). The background increases with respect to the standard selection by a factor 
of two because of the loose cut on Qt < 6MeV/c. Moreover, signal and background 
increase because 

• no cut on the multiplicity two in the SFD search window is needed (gain > 25%). 

• no detector response of the SFD is needed (gain due to single track efficiency of 
95% for two tracks in two planes > 23%.). 

The results are summarized in Tabled The column "Full" shows the results 
obtained with full tracking (cf. section EJ) and background from Monte Carlo. The 
column "MC" shows the result from reconstruction using Monte Carlo backgrounds 
and no vertex fit [HI]. The column "ACC" shows the result of reconstruction using the 
measured accidentals and vertex fit. The column "MC-3 hits" is analogous to column 
"MC" but allowing for 3 hits in the SFD planes [HI] instead of two (cf. section HJ). 
The column "Down" shows the result of using only the downstream detectors and the 
background reconstruction based on accidentals. The integral (0 to 15MeV/c) number 
of events retained after the specific selection cuts of each analysis methods is also given. 
The integral number of measured accidentals (needed for background reconstruction for 
col. "ACC") is about the same as the number of prompt events. The last row of Tabled 
shows the signal fraction extracted from the selected data sample. 
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Table 1. Reconstructed atomic pair events for Q and Ql using full tracking and 
Monte Carlo background (Full), no vertex fit and Monte Carlo background (MC), 
vertex fit and background reconstruction from accidentals (ACC), as "MC" but using 
up to 3 hits instead of two in the SFD acceptance (MC-3 hits), and reconstruction 
only with downstream detectors (Down). Full selected data sample used for signal 
extraction (integral from to 15MeV/c) for each method (selected sample), and ratio 
of extracted signal and selected data sample ( , sl ^ n f ® L , ). 

° r- \ selected sample ' 

Full MC ACC MC-3 hits Down 

Q 5526 ±385 6520 ±370 4670 ±930 8230±440 

Q L 5322 ± 350 6510 ±330 5730 ± 580 8050±380 9280 ± 970 

selected sample - 429000 437060 385870 624880 1503700 



signal Q L 
selected sample 



1.24% 1.49% 1.48% 1.29% 0.6% 



In the case of full tracking (col. "Full" in TableQ) the integrated luminosity was 
(for the Ni-runs in 2001) only 79% of the preselected data due to the availability of the 
MSGC detector. Moreover, these data were not subject to a cut on 2 hits in the SFD 
(section EJ) and thus should be compared with column "MC-3 hits" of Tabled 

First we observe that the number of events in the full sample depends on the 
details of the cut procedures and selection methods. Using downstream detectors only 
(col." down" of Table^) results in a background much larger than the signal increase. 

The consistency of the signal is satisfactory as can be seen from Tabled The signal 
fractions of the selected data sample are the same for the Monte Carlo based method 
and the method based on accidentals (Tabled columns "MC" and "ACC"), and very 
similar for full tracking and MC-3 hits (TableHJ columns "Full" and "MC-3 hits"). 
The large errors of column "ACC" of Tableware due to the limited statistics of the 
measured accidentals. The difference in signal strength for Q and Q L is due to slightly 
different fit regions for Q (Q T < 4MeV/c, \Q L \ < 22MeV/c, Q > 4MeV/c) and Q L 
{Q T < 4MeV/c, 2MeV/c < \Q L \ < 22MeV/c) and for signal searching (Q < 4MeV/c 
and \Ql\ < 2MeV/c, respectively). Thus, the backgrounds in Q and Q L are not 
identical and fluctuations of the background differences in the prompt data and in 
the measured accidental data lead to the (accidentally large) difference in the signal. 
The same argument holds for the Monte Carlo based method. However, there the 
statistical fluctuations of the background are negligible as compared to the measured 
data. Allowing for three hits instead of two in the SFDs yields about the same number of 
events in the signal as were found without using any upstream detector. This indicates 
that the main source of loss in signal is inefficiency in the upstream detectors. 

7. Conclusion 

For the first time a large statistics sample of 7r + 7r~ pairs from atom break-up has 
been detected. Independent tracking procedures and complementary background 
reconstruction strategies lead to compatible results. The background as obtained with 
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Monte Carlo methods is in excellent agreement with the data. We conclude that 
severe systematic errors may be excluded when extracting the atomic pair signal. The 
statistical accuracy of the signal allows for an estimated statistical error on the lifetime 
of the A 27T atom of about 15% In this paper we have analysed only part of our data. 
In view of the full statistics accumulated by the experiment so far and with systematic 
errors estimated to be smaller than the statistical ones J^TJ the goal of the experiment 
to achieve a 10% accuracy for the lifetime of the A 2n atom is in reach. 
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